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It has been clearly established, both experimentally and theoretically, that the 
various regions of the solar-terrestrial system are strongly coupled, that the coupling 
processes exhibit time delays, and that feedback mechanisms exist [cf. Intriligator , 1984; 
Romick et al., 1987; Liu and Edwards, 1988; Roederer, 1988; Schunk, 1989]. For 
example, changes in the solar wind dynamic pressure and the interplanetary magnetic field 
affect the magnetospheric currents and electric fields, which in turn, affect the ionospheric 
convection pattern, electron density morphology, and ion composition at high latitudes. 
The changes in the ionosphere then affect the thermospheric structure, circulation, and 
temperature on a global scale. The changes in the ionosphere-thermosphere system then act 
to modify the magnetospheric processes. The variations in the ionospheric conductivities 
modify the magnetospheric electric fields and the large-scale current system linking the two 
regions. Additional feedback mechanisms occur in the polar cap via the ‘polar wind’ and in 
the auroral zone via ‘energetic ion outflow,’ and these ionospheric ions are a significant 
source of mass, momentum, and energy for the magnetosphere. However, all of the 
coupling and feedback mechanisms have time delays associated with them, which further 
complicates the situation. 

With the above description in mind, the overall goal of our NASA Theory Program 
was to study the coupling, time delays, and feedback mechanisms between the various 
regions of the solar- terrestrial system in a self-consistent, quantitative manner. To 
accomplish this goal, it will eventually be necessary to have time-dependent macroscopic 
models of the different regions of the solar-terrestrial system and we are continually 
working toward this goal. However, with the funding from this NASA program, we 
concentrated on the near-earth plasma environment, including the ionosphere, the 
plasmasphere, and the polar wind. In this area, we developed unique global models that 
allowed us to study the coupling between the different regions. These results are 
highlighted in the next section. 

Another important aspect of our NASA Theory Program concerned the effect that 
localized ‘structure’ had on the macroscopic flow in the ionosphere, plasmasphere, 
thermosphere, and polar wind. The localized structure can be created by structured 
magnetospheric inputs (i.e., structured plasma convection, particle precipitation or 
Birkeland current patterns) or time variations in these inputs due to storms and substorms. 
Also, some of the plasma flows that we predicted with our macroscopic models could be 
unstable, and another one of our goals was to examine the stability of our predicted flows. 

Because time-dependent, three-dimensional numerical models of the solar-terrestrial 
environment generally require extensive computer resources, they are usually based on 
relatively simple mathematical formulations (i.e., simple MHD or hydrodynamic 
formulations). Therefore, another goal of our NASA Theory Program was to study the 
conditions under which various mathematical formulations can be applied to specific solar- 
terrestrial regions. This could involve a detailed comparison of kinetic, semi-kinetic, and 
hydrodynamic predictions for a given polar wind scenario or it could involve the 
comparison of a small-scale particle-in-cell (PIC) simulation of a plasma expansion event 
with a similar macroscopic expansion event. The different mathematical formulations have 
different strengths and weaknesses and a careful comparison of model predictions for 
similar geophysical situations provides insight into when the various models can be used 
with confidence. 
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RESEARCH HIGHLIGHTS 


As a result of our involvement in the NASA STTP program, we published 135 
scientific papers and we presented 161 talks at both national and international meetings. 
These publications and presentations are listed at the back of this report Because of the 
large number of publications, it is not possible to review all of our work in this final report. 
However, in the paragraphs that follow, we highlight some of the studies we conducted, 
with the emphasis on the work done during the last 3-5 years of the program. 


Global Ionosphere 

With support from the NASA Theory Program, we developed the first numerical 
model of the global ionosphere (papers 64, 107, 109, 112, 126, and 133). The model is a 
time-dependent, three-dimensional, high-resolution, multi-ion ionospheric model that 
covers the altitude range from 90-1000 km. With the model, the density distribution of six 
ion species (NO + , 02 + , N 2 + , 0 + , N + , He + ) and the electron and ion temperatures are 
obtained from a numerical solution of the appropriate continuity, momentum, and energy 
equations. 

In recent work (paper 126), the global model was used to study the asymmetry in 
large-scale ionospheric features between the northern and southern hemispheres. The 
comparisons were done for June and December solstice conditions at solar maximum for 
quiet geomagnetic activity. Ideal conditions and diumally-reproducible ionospheric 
features were established in order to elucidate the intrinsic hemispherical differences that are 
associated with the different displacements between the geomagnetic and geographic poles 
and the different atmospheric conditions. In comparing the ionospheric densities in the 
northern and southern hemispheres for a given season, we found the following: (1) The 
winter hemispheres display the most marked universal time (UT) variations due to the 
displacement between the geomagnetic and geographic poles; (2) For similar winter 
conditions, the electron densities in the northern and southern hemispheres can differ by an 
order of magnitude in places at certain UTs; (3) The winter mid-latitude trough is the 
feature that exhibits the largest northern-southern hemisphere difference due to the different 
dipole tilts; and (4) The ‘winter anomaly’ is present in the northern hemisphere at almost all 
UTs, while it is essentially absent in the southern hemisphere. In the northern hemisphere, 
the anomaly maximizes between 1000 and 1600 UT, with the winter/summer peak electron 
density (NmFi) ratio reaching 1.6. In our model, the anomaly results primarily from the 
adopted atmospheric densities. 

Some of the hemispherical differences described above are shown in Figure 1, 
where contours of N e at 300 km are displayed for one particular geometry of solar 
illumination. In the left panel, 1900 UT June solstice conditions are shown, while the right 
panel corresponds to 0700 UT December solstice conditions plotted with an inverted 
latitude axis. The latitude axis was inverted so that summer is at the top and winter at the 
bottom for both the northern and southern hemispheres. 


Global Polar Wind 

The polar wind is an ambipolar outflow of thermal plasma from the terrestrial 
ionosphere at high latitudes. The outflow typically begins at about 800 km. As the 
ionospheric ions flow up and out of the topside ionosphere along diverging geomagnetic 
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Figure 1. Global plot of the electron density variation at 300 km in a magnetic latitude-MLT 
reference frame. The left panels are for June solstice at 1900 UT and the right, panels are for 
coicHrp at 0700 ITT The too half of the figure is for the summer hemisphere and the 

in the southern hemisphere for a larger offset between die magnetic and geographic poles. The Ne 
contours are labeled in logio Ne (cm*3). From Paper 126. 
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field lines, they are accelerated and eventually become supersonic (above about 1300 km). 
Although the classical polar wind had been studied for twenty years, all of the studies were 
based on either steady state or time-dependent, one-dimensional models applied to a single 
location. However, we constructed the first three-dimensional, time-dependent, multi-ion 
model of the global polar wind in order to study the temporal evolution of ion outflow 
during magnetic storms and substorms (paper 130). The model covers the altitude range 
from 120 to 9000 km. At low altitudes (120-800 km), three-dimensional distributions for 
the NO + , 02 + , N 2 + , N + , and 0 + densities and the ion and electron temperatures are 
obtained from a numerical solution of the appropriate continuity, momentum, and energy 
equations. At high altitudes (500-9000 km), the time-dependent, nonlinear, hydrodynamic 
equations for 0+ and H + are solved self-consistently with the ionospheric equations taking 
into account collisions, charge exchange chemical reactions, flux tube divergence, and ion 
temperature anisotropies. The model can describe supersonic ion outflow, shock 
formation, and ion energization during plasma expansion events. 

In the first application of the model, we studied the temporal response of the global 
polar wind to changing magnetospheric conditions. During the simulation, the magnetic 
activity level changed from quiet to active and then back to quiet again over a 4. 5-hour 
period. When the activity level increased, the auroral oval expanded, particle precipitation 
became more intense, and the convection electric field strengths increased. The reverse 
occurred during declining magnetic activity. From this study we found the following: (1) 
Plasma pressure changes in the ionosphere due to T e , Tj, or electron density variations 
produce perturbations in the polar wind. In particular, plasma flux tube motion through the 
auroral oval and high electric field regions produces transient large-scale ion upflows and 
downflows. At certain times and in certain regions both inside and outside the auroral 
oval, H + -0 + counterstreaming can occur; (2) The density structure in the polar wind is 
considerably more complicated than in the ionosphere because of both horizontal plasma 
convection and changing vertical propagation speeds due to spatially varying ionospheric 
temperatures; (3) During increasing magnetic activity, there is an overall increase in T e , Ti, 
and the electron density in the F region, but there is a time delay in the buildup of the 
electron density that is as long as five hours at high altitudes; (4) During increasing 
magnetic activity, there is an overall increase in the polar wind outflow from the 
ionosphere, while the reverse is true for declining activity. In certain regions, however, 
localized ionospheric holes can develop during increasing magnetic activity, and in these 
locations the polar wind outflow rate is reduced; and (5) During changing magnetic activity, 
the temporal evolution of the ion density morphology at high altitudes (~ 9000 km) can be 
opposite to that at low altitudes (~ 500 fan). 

Some of these results are seen in Figure 2, which shows ‘snapshots’ of the H + 
density (logiofl) at 500 km (bottom row), 2500 km (middle row) and 9000 km (top row) 
above the polar region for times corresponding to increasing activity (fi), high activity (( 4 ), 
and declining activity fo). At the start of increasing magnetic activity (/ 1 ), the H + density 
morphology is similar to the 0 + morphology (not shown), with relatively high densities in 
the oval in the midnight-dawn-noon sector, moderate densities in the central polar cap, and 
low densities in the dayside trough region. This result is not surprising, since at low 
altitudes 0 + is the source of H + via the 0 + + H <=> H + + O reaction. As the storm 
intensifies O4), there is a large-scale H + upwelling that occurs in association with both the 
electron heating in the auroral oval due to increased particle precipitation and the ion heating 
in the regions of large electric fields. During declining magnetic activity the H + 
density at 500 km is generally higher than what it was at any time during the storm. 'Hiis 
time delay in the H + density buildup at low altitudes is simply a consequence of the time 
delay in the 0 + density buildup (not shown). However, at 2500 km, there are some polar 
regions where the H + density is lower than at the storm’s main phase (fy) and other polar 
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Figure 2. H+ density contours at 500 km (bottom row), 2500 km (middle row) and 9000 km (top 
row) for times corresponding to increasing (?i), high (f4) and declining (t6) magnetic activity. A 
magnetic latitude-MLT reference frame is used. Note that different color scales are associated with 
the different altitudes. From paper 130. 
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regions where it is higher. Furthermore, at high altitudes (9000 km), there is an overall 
reduction in the H + density at t(, relative to the storm’s main phase fo) even though there is 
an overall H + density increase at low altitudes (500 km). The reason for the difference in 
the H + density morphologies at high and low altitudes during changing magnetospheric 
conditions is that there is a finite propagation time for low-altitude perturbations to reach 
high altitudes. The situation is further complicated by horizontal plasma convection, 
changing vertical propagation speeds with changing plasma temperatures, and localized 
regions where H + downflows can occur. 


Global Plasmasphere 

During the past two decades there have been several theoretical investigations of the 
coupling between thermal plasmas in the ionosphere and plasmasphere [e.g., Horwitz, 
1983; Lemaire, 1985; and references therein]. Yet, there remains many unresolved 
questions in fundamental areas of research relating to the plasmasphere. An example of 
one of these areas is die location of the plasmapause. In spite of the obvious importance of 
the plasmapause, there remains an uncertainty as to the physical processes involved in its 
formation and location. For many years it has been accepted that the location of the 
plasmapause is controlled by a dynamic process associated with E x B convection in the 
magnetosphere [e.g., Nishida , 1966; Brice , 1967]. However, Lemaire [1985] has listed 
several problems with this model. 

In order to study ionosphere-plasmasphere coupling phenomena and plasmapause 
formation mechanisms, we developed a three-dimensional, time-dependent model of the 
global plasmasphere. The model solves the nonlinear hydrodynamic continuity and 
momentum equations along closed flux tubes for multiple H + streams arising from the 
conjugate hemispheres. The model takes into account the production and loss of H + , 
collisions of H + with 0 + , and the motion of plasmaspheric flux tubes in response to 
convection electric fields (cross-L drifts). The model retains the inertial term in the 
momentum equation so that supersonic flow and various low frequency wave phenomena 
can be studied. 

In our first study of the global plasmasphere (paper 132), the model was run for 
solar minimum conditions until diumally reproducible results were obtained, indicating that 
the plasmasphere was fully filled. Equatorial densities obtained during this study are 
shown in the top panel of Figure 3. An interesting feature seen is the diurnal variation in 
density. Earthward of approximately L - 2.5, a minimum in density occurs at local 
midnight due to decreasing ionization levels in the ionosphere. However, the density of 
plasma near the plasmapause maximizes at local midnight due to a compression of the 
plasma as flux tubes convect toward the Earth in this time sector. The effect of expansion 
of the plasma can be seen in the dusk region, where the density of plasma converting near 
the plasmapause is a minimum. Another view of the global plasmasphere is seen in the 
bottom panel of Figure 3, where H + densities are shown in the dusk-dawn meridian plane. 
This panel clearly shows the dusk bulge, indicated by the increased volume of the 
plasmasphere in the dusk sector relative to the dawn sector. When the results of the global 
model were compared with applicable measurements, it was found that the model results 
were within a factor of two of both the whistler and satellite observations. Near the 
plasmapause, changes in flux-tube volume due to cross-L drifts led to a factor of 3 
variation in equatorial density, with the highest densities occurring near local midnight 
where the flux-tube volume was lowest 
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Figure 3. Contours of H + density in the equatorial plane of the magnetosphere (top). L value 
is indicated on the radial axis and local time on the outer circle. Sunlit and dark portions of the 
Earth are shown in the center of the figure. In the bottom panel, contours of density in the 
plasmasphere are shown at two different times. Black lines are plotted at Zr-shells of 1.5 and 3.0 and 
the outer red line is plotted at L — 5.0. The Earth is shown in light blue in the center. For both 
panels, the color scale is logj 0 nfH"* - ) in cm -3 and portions of the figure shown in black represent 
the plasma trough region outside of the plasmasphere. From Paper 132. 



9 


Polar Electrodynamics 

In studies involving numerical simulations of the ionosphere, polar wind, 
plasmasphere, and thermosphere, the magnetospheric processes (field-aligned currents, 
convection electric fields, and particle precipitation) are usually taken as known inputs [cf. 
Roble etal., 1988; Fuller -Rowell and Rees, 1983; papers 109, 130, 132, and 133]. These 
inputs are obtained from either empirical models or direct measurements. However, the 
field-aligned currents that flow into and out of the ionosphere are connected via horizontal 
currents that flow in the lower ionosphere (E region). These large-scale currents, the 
auroral conductivity enhancements due to precipitating energetic electrons (upward 
currents), and the convection electric fields are not independent, but instead, are related via 
Ohm’s law and the current continuity equation. In this regard, we obtained a numerical 
solution of the appropriate equations in older to study the current-conductivity-electric field 
coupling between the ionosphere and magnetosphere (papers 95, 111, and 116). Our 
polar electrodynamics model is similar to the models developed earlier by others [cf. 
Kamide et al., 1986; Ahn et ai, 1986; Richmond and Kamide, 1988]. With a measured 
field-aligned current distribution and conductivities calculated from our global ionospheric 
model, we are able to calculate the self-consistent convection electric fields and horizontal E 
region currents. We are also able to take into account the F region contribution to the 
ionospheric conductivity. 

The current-conductivity-electric field coupling at high latitudes is strongly 
dependent on the direction of the interplanetary magnetic field (IMF). When the IMF is 
southward, the polar region is characterized by a two-cell pattern of plasma convection, 
Region 1 and 2 field-aligned current systems, and enhanced conductivities in the auroral 
oval. The overall consistency of these features is shown in Figure 4, where ionospheric 
conductivities, electric potentials, field-aligned currents, and E region currents are shown 
for one of our model simulations with B z < 0 and By <0 (paper 1 16). For northward 
IMF, on the other hand, the situation is more complicated and our work on this subject will 
be discussed later. 


Large-Scale Ionospheric Structures 

During the last decade, significant progress has been made in understanding the 
formation and evolution of the major morphological features observed in the high latitude 
ionosphere. However, the ionosphere exhibits a considerable amount of spatial structure 
that is not reproduced with the global numerical models when ‘fixed’ empirical models are 
adopted for the magnetospheric convection and precipitation patterns. There are small-scale 
(5 1 km), medium-scale (~ 10 km), and large-scale (> 100 km) density structures [cf. 
Buchau et al., 1985; Weber et al., 1986; Robinson et al., 1985; Tsunoda, 1988]. The 
small-scale structures appear to be produced within and on the edges of the larger structures 
through plasma instabilities and are typically referred to as density irregularities [Kelley et 
al., 1982; Lee, 1984]. The medium and large-scale density structures have been called 
‘blobs’, ‘patches’, and ‘enhancements’. Relative to background densities, the 
perturbations associated with the medium and large-scale structures vary from about 10% 
to a factor of 100. 

With regard to the medium and large-scale density structures, numerous structuring 
mechanisms have been proposed, including structured precipitation, structured electric 
fields, structured heat flow, transport from distant sources, structured enhanced decay, as 
well as time varying processes of all kinds [cf. Kelley et al., 1982; Heelis et al., 1986; de la 
Beaujardiere et al., 1985; Weber et al., 1984; Robinson et al., 1985; paper 119]. Although 
all of the proposed mechanisms can produce density structures in a qualitative sense, very 
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Figure 4. Pedersen conductivity (top left panel), Birkeland currents (top right panel), ionospheric 
electric potential (bottom left panel) and ionospheric currents (bottom right panel) for B x < 0 
and B v < 0. The difference between contour levels is 0.1 mho for the conductivity plot. For 
the Birkeland currents, current flowing into the ionosphere is denoted by solid lines, while current 
flowing away from the ionosphere is shown by dashed lines (the contour separation is 0.2 /i A/m^). 
For the electric potential, solid lines denote positive potentials and dashed lines denote negative 
potentials (the contour separation is 6 kV). For the horizontal current, the vector scale is 400 A/km. 
From Paper 116. 
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little attention has been given to detailed quantitative studies of structure formation 
mechanisms. However, it is important to understand the production and transport 
mechanisms of ionospheric density structures because plasma enhancements affect the 
momentum and energy coupling between the ionosphere and thermosphere. Consequently, 
we conducted a systematic quantitative study of the origins, lifetimes, and transport 
characteristics of medium and large-scale density structures using our time-dependent, 
three-dimensional ionospheric model. We studied the formation of plasma ‘blobs’ or 
‘patches’ due to particle precipitation in the dayside oval (papers 72 and 88), the effect that 
structured electric fields have on the ionosphere (paper 103), the extent to which large 
magnetic-field-aligned ion drifts affect density structures (paper 104), the effect of time- 
varying convection and precipitation patterns (papers 35 and 52), and the lifetimes and 
transport characteristics of density structures for different seasonal, solar cycle, and IMF 
conditions (papers 101 and 119). 

The above studies were particularly useful for elucidating basic structuring 
processes. In particular, we found that: (1) Time- varying plasma convection and particle 
precipitation patterns can readily produce ionospheric structure on a variety of scales; (2) 
The observed electron precipitation in the polar cap can account for ‘blob’ creation if the 
plasma is exposed to the precipitation for 5-10 minutes; (3) Structured electric fields 
produce structured electron densities, ion temperatures, and ion composition; (4) The 
lifetime of an F region density structure depends on several factors, including the initial 
location where it was formed, the magnitude of the perturbation, season, solar cycle and 
IMF; and (5) Depending on the IMF, horizontal plasma convection can cause an initial 
structure to break up into multiple structures of various sizes, remain as a single distorted 
structure, or become stretched into elongated segments. 

In order to model a convection pattern with electric field structure superimposed, 
measurements are needed not only along the satellite tract, but also in a direction 
perpendicular to the track, so that the entire spatial extent of the structure can be 
determined. Unfortunately, the required data set was not available and the model study 
(paper 103) had to be based on a semi-empirical description of electric field structure. The 
structure was assumed to be associated with a discrete current system into and out of the 
polar ionosphere; the current sheets for each structure were assumed to be sun-aligned and 
had a finite length. The potential disturbance associated with the structured current sheets 
was assumed to be a Volland two-cell pattern [ Volland , 1978], as shown in Figure 5a. 
Each electric field structure was described by one such 2-cell pattern, but the spatial size of 
the structure and potential drop across it varied from structure to structure. In the model 
study, a total of 45 structures were randomly distributed in the polar cap. 

The structured potential distribution was added to a standard, large-scale Volland 
potential distribution, which was used to simulate the background convection pattern. In 
the model study, eight representative convection trajectories were selected and flux tubes of 
plasma were followed as they convected along these trajectories across the polar cap in an 
anti-sunward direction. The calculations were done both with and without the electric field 
structure. The top panel of Figure 5b shows the trajectories for the unperturbed or 
background convection pattern, while the bottom panel shows the resulting trajectories 
when the structured electric fields are superimposed on the background pattern. The 
electric fields that would be measured by a satellite crossing these convection patterns along 
the dusk-dawn magnetic meridian are shown in Figure 5c. The top panel shows the electric 
field signature of a classical two-cell convection pattern, while the bottom panel shows the 
kind of structure that is frequently observed when the IMF is northward. 

In the calculations, the trajectories were all started from daytime steady state 
conditions at the locations marked by the dot in Figure 5b, and the response of the plasma 
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Figure 5a. Model of an electric field structure. 
Electric potential drawn at 1 kV intervals are shown 
in a magnetic x-y coordinate system. The potential 
varies from —6 to 6 kV for this structure. 
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Figure 56. Plasma trajectories for a 2-cell un- 
perturbed convection pattern (top panel) and for 
a structured convection pattern (bottom panel) in 
the MLT-magnetic latitude reference frame. Circles 
denote the start location for each trajectory, while 
tick marks indicate one-hour time intervals. 
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Figure 5c. Electric field variation along the dusk- 
dawn magnetic meridian for the unperturbed 2-cell 
convection pattern (top panel) and for the struc- 
tured convection pattern (lower panel). The elec- 
tric field is given in mV/m. 
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Figure 5 d. N m F 2 variation for each of the trajecto- 
ries shown in Figure 56 plotted versus the flux tubes 
x location. The unperturbed (top) and perturbed 
(bottom) cases are shown. log 10 N m F 2 (cm" 3 ) is 
plotted. From Paper 103. 
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was calculated as the flux tubes followed their convection paths. Figure 5d shows how 
N m F 2 varies along each trajectory for both the unperturbed (top panel) and structured 
(bottom panel) convection patterns. Note that the densities are plotted versus the trajectory 
jt-axis (see Figure 5a). For the unperturbed pattern, the plasma following the eight 
trajectories convects across the polar cap and eventually into darkness, during which time 
N m F 2 decays. However, the important feature to note is that the N m F 2 variations for all 
eight trajectories are very similar. In sharp contrast, the N m F 2 variations along the eight 
trajectories that contain structured electric fields are very different. Upon reaching the 
nightside auroral zone, the NmF2 values for the different flux tubes can differ by as much 
as an order of magnitude. The differences are a result of different transpolar crossing times 
and different ‘localized’ electric field strengths. Localized enhancements in the electric field 
produce ion temperature enhancements and increased 0+ — » N0 + conversion rates [Schunk 
etal., 1975], and this, in turn, yields lower N m / r 2 values. 


Northward IMF Features 

It is well known that the electric fields, particle precipitation, auroral conductivity 
enhancements, and Birkeland currents that couple the magnetosphere-ionosphere system 
are strongly dependent upon the direction of the IMF. When the IMF is southward, the 
Birkeland (or field-aligned) currents flow in the Region 1 and 2 current sheets, the F region 
plasma convection exhibits a two-cell structure with anti-sunward flow over the polar cap, 
the conductivity enhancements are confined to the statistical auroral oval, and the auroral 
electron precipitation is also confined to the classical oval. However, when the IMF is 
northward, the situation is considerably more complicated and less clear. In this case, an 
additional field-aligned current system occurs in the polar cap called the NBZ currents 
[Maezawa, 1976; Iijima et al., 1984]; plasma convection can be sunward in the polar cap 
and the pattern can assume multi-cell, severely distorted two-cell or turbulent characteristics 
[Burke et al., 1979; Potemra et al ., 1984; Frank et al., 1986; Heppner and Maynard, 
1987]; and particle precipitation occurs in the polar cap which can be uniform, in the form 
of multiple sun-aligned arcs, or in a 0-aurora configuration [Buchau et al., 1983; Frank et 
al., 1986]. 

We conducted a number of theoretical studies in an effort to elucidate the 
ionospheric conditions for northward IMF. Using an E region conductivity model we 
developed (paper 111) and our polar electrodynamic model, we conducted the first current- 
conductivity simulations for northward IMF. In our first paper, we considered the 
magnetosphere to be both a ‘constant current source’ and a ‘constant voltage source’ and 
determined the conditions leading to idealized two-cell, three-cell, and four-cell plasma 
convection patterns (paper 95). When the various cases were compared with 
measurements, our results indicated that the magnetosphere acts more like a voltage source 
than a current source. 

Subsequently, Iijima and Shibaji [1987] presented empirical models of Birkeland 
currents for different IMF By and B z conditions, and we used these current patterns in our 
electrodynamic simulations in an effort to determine whether three-cell, four-cell, or 
distorted two-cell convection patterns occur for northward B z (paper 1 16). Unfortunately, 
because of a data gap on the dayside at latitudes equatorward of about 80° magnetic, the 
measured NBZ currents had to be extrapolated in the noon sector. Depending on how the 
currents were extrapolated, different convection patterns could be obtained. Figure 6 
illustrates the above problem for the case when the IMF is strongly northward (B z > By > 
0). Shown in the upper-left dial of this figure are the unmodified Birkeland currents 
measured by Iijima and Shibaji [1987]. The upper-right dial, on the other hand, shows a 
possible extension of these currents into the dayside where measurements are lacking. The 
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Figure 6. Birkeland currents (upper panel), ionospheric electric potential (middle panel) and 
ionospheric currents (lower panel) for B z > B v > 0. The upper-left dial represents the unmodified 
Birkeland currents of Iijima and Shibaji [1987], while the currents in the upper-right dial have been 
extended into the dayside and nightside. In the upper panel, current flowing into the ionosphere is 
plotted as solid curves, while current flowing away from the ionosphere is plotted as dashed curves 
(the difference between contour levels is 0.3 fiA/m?). In the middle panel the difference between 
contour levels is 3 kV. For the lower panel the vector scale is 400 A/km. From Paper 116. 
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extension was accomplished by interpolating along lines of constant latitude so that the 
dusk Region 1 system connects to the dawnside NBZ system. The electric potential 
patterns calculated from the two Birkeland current patterns are shown in the middle panel. 
In comparing the two potential patterns, it is apparent that the extension of the current on 
the dayside is important for northward IMF. For the unmodified currents on the left-hand 
side, the flow in the morning sector is broken into multiple cells. Although there is 
evidence for three-cell convection, the NBZ currents are so strong in comparison to the net 
Region 1/Region 2 currents, that the strongest flow occurs within the polar cap and the 
rotation of the equatorward cell in the dawn sector is reversed from its normal sense. For 
the dayside extension of the currents in the upper-right dial, the flow also has a three-cell 
nature, but because of the connection of flow on the dayside, convection appears more as a 
distorted two-cell pattern similar to that shown by Heppner and Maynard [1987]. Although 
the convection pattern described above supports the distorted two-cell configuration 
suggested by Heppner and Maynard [1987], different extrapolations of the NBZ current 
system in the noon sector can produce both three-cell and four-cell convection patterns. 

Once the convection patterns for northward IMF are identified, we have the 
additional capability of determining their effect on the F region. For example, in one study 
the idealized multi-cell plasma convection patterns described above for northward IMF 
were used in conjunction with our time-dependent, three-dimensional ionospheric model in 
order to study the characteristic F region signature associated with multi-cell convection 
patterns (paper 86). As expected, there are major distinguishing F region features 
associated with the different convection patterns, particularly in the polar cap in winter 
(Figure 7). For two-cell convection, the anti-sunward flow of plasma from the dayside 
into the polar cap acts to maintain the densities in this region. For four-cell convection, on 
the other hand, the two additional convection cells in the polar cap are in darkness most of 
the time, and the resulting 0 + density decay acts to produce twin polar holes that are 
separated by a sun-aligned ridge of enhanced ionization due to 0-aurora precipitation. For 
three-cell convection, only one polar hole forms (not shown). 


Stability of Macroscopic Flows 

Numerous mathematical formulations have been used over the years to describe 
plasma flows in the solar-terrestrial environment, including magnetohydrodynamic, kinetic, 
semi-kinetic, hydromagnetic, generalized moment, and hydrodynamic formulations. All of 
these formulations have intrinsic limitations when applied to three-dimensional macroscopic 
flows. For example, the kinetic and semi-kinetic formulations are usually applied to steady 
state situations, and hence, wave phenomena are not included. Also, although time- 
dependent MHD and hydrodynamic models allow for wave excitation, only low frequency 
waves are possible. Hence, in general, the microphysics (non-Maxwellian distributions, 
plasma instabilities, wave-particle interactions, etc.) is not included in full-scale global 
simulations. 

Unstable flow conditions can occur in a variety of situations. For example, the 
edges of convecting F region blobs and patches are susceptible to gradient-drift instabilities 
[cf. Tsunoda, 1988], but these instabilities generally don’t affect the large-scale flow. 
Also, in regions of rapid E x B convection, die ion velocity distributions in the E and F 
regions can become highly non-Maxwellian [cf. St.-Maurice and Schunk, 1979; Hubert, 
1983; Lockwood et al, 1987] and the resulting instabilities could affect ionosphere- 
thermosphere momentum and energy coupling. Likewise, the large-scale horizontal 
currents that flow in the E region (see Figure 4) are susceptible to Farley-Buneman 
instabilities [cf. St.-Maurice , 1985; Haldoupis, 1989]. Finally, we note that the polar wind 
distributions predicted by the macroscopic models are highly non-Maxwellian and 
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potentially unstable, particularly when energetic beams and conics pass through the polar 
wind. Therefore, when macroscopic models are used to predict plasma flow conditions, it 
is important to study the stability of the predicted flows. Part of our NASA Theory 
Program was devoted to such an effort and in the following paragraphs we highlight our 
polar wind instability studies. 

The classical polar wind corresponds to an outflow of thermal plasma from the high 
latitude ionosphere. Above about 3000 km, the H + flow becomes supersonic and 
collisionless, and the H + velocity distribution function becomes highly non-Maxwellian. 
The non-Maxwellian features include a temperature anisotropy with T/ 11 > Tj 1 and an 
asymmetry in the form of an elongated tail in the upward direction. At a distance of 10 
earth radii, Tj 11 / T^- ~ 50. Since this highly non-Maxwellian H + velocity distribution is 
potentially unstable and since the resulting wave-particle interactions could have a 
significant effect on the flow of mass, momentum, and energy between the polar 
ionosphere and magnetosphere, we studied the linear stability of the polair wind with regard 
to the excitation of electrostatic waves (paper 93). Despite the highly non-Maxwellian 
character of the H + distribution function, we found that the ‘classical’ polar wind is 
remarkably stable for a wide range of conditions. 

Although the classical polar wind occurs the bulk of the time, during increasing 
magnetic activity and for certain flux tubes that pass through the cusp, die classical picture 
of the polar wind may not be appropriate. Figure 8 shows schematically a possible 
mechanism by which the polar wind deviates from the classical picture due to plasma 
convection. As the flux tube convects into the cusp region, the ions at low altitudes 
undergo parallel and perpendicular energization. These energetic ions eventually catch up 
with the slower “classical” polar wind at high altitudes as the flux tube convects into the 
polar cap region. The stability of a cusp-generated H + beam passing through the classical 
polar wind was studied in paper-123. The stability of the plasma was studied for a range of 
electron-ion temperature ratios (0.1 <, T e (Ti £ 10) and beam-to-background ion density 
ratios (0.1 <> ntJint, + /»/) £ 0.9). Table 1 summarizes the results found in our study. For 
most of the cases, the polar wind is unstable for fairly small ion beam-plasma relative 
drifts. The plasma is more unstable for larger electron temperatures and for comparable ion 
beam and background ion densities (rib ~ «»). For large electron temperatures (T e = 10 7/), 
the ion/ion plasma-acoustic instability is the first mode to be triggered, and the waves that 
propagate parallel to the magnetic field are destabilized first. For lower electron 
temperatures, the ion/ion cyclotron modes tend to be triggered first 


Auroral Dynamics 

We used our 2.5-dimensional particle-in-cell (PIC) code to study auroral plasma 
processes (papers 90 and 91). In particular, we conducted a systematic study of the plasma 
processes driven by a field-aligned current sheet with a finite thickness perpendicular to the 
B field. Both wide and narrow (relative to the ion gyroradius) current sheets were 
considered. In these simulations, we observed double layers, upflowing ions in the form 
of beams and conics, both upflowing and downflowing accelerated electron beams, several 
wave modes (EIC, lower hybrid, VLF, ELF, and HF), and associated nonlinear plasma 
phenomena. 
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Figure 8. Schematic diagram showing how ion acceleration in the cusp coupled with horizontal 
plasma convection can lead to an energetic ion population passing through the classical polar wind. 
From Paper 123. 


Table 1. H + beam-plasma normalized critical drift velocity- 
above which the plasma is unstable for different electron-to- 
ion temperature ratios and different beam-to-background ion 
density ratios for Ti/Ti = 0.1. The critical drift velocity is 
normalized with respect to the parallel thermal speed of the 
background ions. From Paper 123. 
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n b /(n b + n i) 


0.1 

0.5 

0.9 
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> 11 

3.4 

> 11 

1.0 

2.7 

0.71 

1.30 

10 

0.3 

0.18 

0.26 
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Subsequently, our attention was focused on studying the effect that auroral 
processes have on the ionospheric densities and temperatures. In paper 97, we studied the 
effect of ionospheric return currents on auroral electron temperatures. This study was 
motivated by the fact that in a partially-ionized plasma an electron heat flow can occur in 
response to either an electron temperature gradient (thermal conduction) or an electron 
current (thermoelectric heat flow). The former process has been extensively studied, while 
the latter process has received relatively little attention. Therefore, we used our time- 
dependent, three-dimensional model of the high-latitude ionosphere to study the effect of 
field-aligned ionospheric return currents on auroral electron temperatures for different 
seasonal and solar cycle conditions as well as for different heat fluxes at our upper 
boundary (800 km). 

Figure 9 shows electron temperature profiles for three values of the field-aligned 
return current for summer and winter conditions at solar minimum. The field-aligned 
currents are 0 (solid curves), —1 x 10'^ (dotted curves), and —5 x 10'^ (dashed curves) amp 
nr 2 . The profiles were calculated with the total heat flow through the upper boundary set 
to zero. The solid curves correspond to the typical case of no field-aligned current and no 
heat flow through the upper boundary. Consequently, the electron temperature is constant 
at high altitudes owing to the dominance of thermal conduction. For all cases shown in 
Figure 9, thermoelectric heat transport corresponds to an upward flow of energy at all 
altitudes, but it can be a source or sink of heat depending on the shape of the electron 
density profile. For example, if we consider the solar minimum- winter case with /n = -5 x 
10‘ 5 amp nr 2 , thermoelectric heat transport is a sink below 200 km, a source between 200 
and 280 km, and a sink above 280 km. The effect of an increased magnitude of the return 
current is merely to enhance this result 


Plasma Expansion Phenomena - 

We also studied plasma expansion phenomena. Based on the pioneering work of 
Gurevich et al. [1966] and subsequent work [cf. Samir et al., 1983], it is now well known 
that during a rapid plasma expansion, the ion distribution functions can depart significantly 
from a Maxwellian, and consequently, superthermal tails, plasma instabilities and wave- 
particle interactions may be important. From the macroscopic modelling viewpoint it is 
important to understand the basic physics involved in a plasma expulsion so that reliable 
macroscopic models can be developed for a range of solar-terrestrial applications (solar 
wind, polar wind, interhemispheric flow, etc.). 

Our work in this area involved both small-scale and macroscopic plasma 
expansions. Initially, we obtained a numerical solution of the time-dependent Vlasov- 
Poisson system and studied single-ion (H + ), multi-ion (H + , 0 + ), and counterstreaming 
(H + -H + ) plasma expansions (papers 26, 33, and 42). These small-scale simulations were 
for collisionless , one-dimensional plasmas, and hence, were relevant to plasma expansions 
along a magnetic field. However, since the B field was not included in these simulations, 
the excitation of certain wave modes (EIC, lower hybrid, etc.) was excluded. 
Consequently, we then conducted 2.5-dimensional PIC simulations for collisionless 
plasma expansions along a B field so that the results could be compared to the 1-D cases 
(paper 79). We also conducted macroscopic (hydrodynamic) plasma expansions to model 
conditions in the collisionless polar wind (papers 68, 81, and 115). Again, these studies 
were for expansions along the B field. We found that many of the expansion 
characteristics seen in the small-scale simulations were also seen in the macroscopic 
simulations, but the macroscopic simulations produced additional results not contained in 
the small-scale simulations and vice versa. 
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Figure 9. Electron temperature profiles in the nocturnal auroral oval for three values of the field- 
aligned return current at solar minimum for both summer and winter conditions. The field-aligned 
current values are 0 (solid curves), —1 x 5 ^ (dotted curves), and —5 X 10 ® (dashed curves) amp 
m“ 2 . From Paper 97. 
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Our most recent work involved macroscopic plasma expansions across a magnetic 
field including the effect of collisions (paper 134). TTie model we developed is based on a 
numerical solution of the time-dependent, nonlinear, coupled continuity and momentum 
equations for the expanding ions, background ions, and electrons. In this model, the 
equations are height-integrated so that only plasma expansions across the geomagnetic field 
can be modelled (i.e., 2-D model). In a model simulation, a neutral gas or plasma is 
injected into the ionosphere and the subsequent temporal evolution is followed. The model 
takes into account ionization processes, chemical reactions, pressure gradients, ion inertia, 
stress tensor (collisional viscosity and finite Larmor radius effect), gravity. Coulomb 
collisions between the different charged particles, ion collisions with both the background 
and injected neutrals, and polarization electric fields in both the expanding and background 
plasmas. Figure 10 shows an example of the 2-D simulation results. This figure compares 
the expansion characteristics for gases released at 200 and 400 km. The two cases are for 
similar conditions; 10 kilogram Barium release with an initial cloud radius of 2 km and a 
release velocity of 5 km/s perpendicular to B, which is directed out of the figure. Shown 
in Figure 10 are contours of the ion cloud Pedersen conductivity and the ion velocity field 
at t = 1, 10, and 25 seconds after the gas release. Also shown in the figure is the motion of 
the neutral Ba cloud, which can be clearly seen at 400 km. The neutral cloud is shown by a 
single circle that expands as it moves in the released gas direction (upward in the figure). 
In comparing similar Ba releases at 400 and 200 km, the following features are evident: (1) 
at both altitudes, the ion clouds skid across magnetic field lines; (2) the expansion of both 
the ion and neutral clouds is slower at low altitudes than at high altitudes; (3) the cloud 
Pedersen conductivity is about 200 times greater at 200 km than at 400 km; and (4) the ion 
cloud expansion is more symmetric at 200 km than at 400 km. 


Solar Wind/Interplanetary Medium 

Recently, we used generalized transport equations for both subsonic flow along 
SAR-arc field lines (paper 94) and for the supersonic polar wind (papers 98 and 120). 
These bi-Maxwellian-based 16-moment transport equations have an advantage over both 
the standard collision-dominated transport equation (Chapman and Enskog equations) and 
the collisionless kinetic equations in that they can describe both collision-dominated and 
collisionless regimes and provide a continuous transition between them. In the collision- 
dominated regime, they reduce to the usual diffusion and thermal conduction equations, 
while in the collisionless regime they describe temperature anisotropies and collisionless 
heat flow to the same level as kinetic equations. Another advantage is that velocity 
distribution functions can be readily obtained at any point in space. As a spinoff of this 
work, we applied the equations to the solar wind and were able to obtain solutions from 30 
R s (solar radii) to the orbit of the Earth at 215 R s (papers 131 and 135). We found that the 
bi-Maxwellian based 16-moment equations can properly describe the evolution of the 
proton distribution function [Marsch etal., 1982]. 

In order to achieve our ultimate goal of obtaining a numerical model of the solar- 
terrestrial system from the solar corona to the lower thermosphere, we developed both 2-D 
and 3-D hydrodynamic models and a 2-D MHD model of steady corotating solar wind 
streams (papers 1 17 and 118). All three models were used to study the latitudinal structure 
of the solar wind. We found that with the MHD model the latitudinal structure of the 
steady solar wind is decided by a dynamic balance between the magnetic and gas pressures. 
Our calculations also indicated the existence of a proton number density maximum around 
the magnetic neutral line regardless of whether there was a density maximum or minimum 
at our inner boundary at 0.15 AU. 




Figure 10. Contours of the ion cloud Pedersen conductivity and the ion velocity field at t = 1, 10, 
and 25 seconds after Barium releases at 400 km (top panel) and 200 km (bottom panel). The B 
field is directed out of the figure. The contour separation is 10~ 3 mho in the upper panel and 0.5 
mho in the lower panel. From Paper 134. 
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47. R.W. Schunk, High latitude ionospheric predictions, presented at the GISMOS 

Workshop, October 24-26, 1984; Boulder, Colorado. 

48. R. W. Schunk, Global models of ionospheric parameters, presented at the GISMOS 

Workshop, October 24-26, 1984; Boulder, Colorado. 

49. R. W. Schunk, Ionosphere-magnetosphere coupling. Invited Talk, presented at the 

Solar-Terrestrial Theory Program Meeting/Workshop, November 14-16, 1984; 
Los Alamos, New Mexico. 

50. R. W. Schunk and J. J. Sojka, Critical tests of a global ionospheric model. Invited 

Talk, AGU Fall Meeting, San Francisco, California; EOS, 65, 1026, 1984. 

51. R.W. Schunk, The flow of plasma in the solar-terrestrial system. Colloquium, 

Department of Physics and Astronomy, Brigham Young University, January 9, 
1985; Provo, Utah. 

52. N. Singh, H. Thiemann, and R. W. Schunk, Electric fields in 2-D double layer 

simulations and their comparison with satellite observations, presented at the 2nd 
International School for Space Simulation, February 4-15, 1985; Kauai, Hawaii. 

53. N. Singh and R. W. Schunk, Temporal behavior of density perturbations in the polar 

wind, presented at the 2nd International School for Space Simulation, February 
4-15, 1985; Kauai, Hawaii. 

54. C. E. Rasmussen, R. W. Schunk, and J. J. Sojka, Modeling of the polar ionosphere, 

presented at the 2nd International School for Space Simulation, February 4-15, 
1985; Kauai, Hawaii. 

55. A. R. Barakat and R. W. Schunk, Effect of hot electrons on the polar wind, 

presented at the 2nd International School for Space Simulation, February 4—15, 
1985; Kauai, Hawaii. 

56. J. J. Sojka and R. W. Schunk, A theoretical study of the global F region for June 

solstice and low magnetic activity, presented at the 2nd International School for 
Space Simulation, February 4—15, 1985; Kauai, Hawaii. 

57. N. Singh, H. Thiemann, and R. W. Schunk, Ion conics as seen in simulations of 2- 

D potential structures, presented at the Chapman Conference on Ion Acceleration 
in the Magnetosphere and Ionosphere, June 3-7, 1985; Wellesley, 
Massachusetts. 

58. N. Singh and R. W. Schunk, Ion acceleration in expanding ionospheric plasmas, 

presented at the Chapman Conference on Ion Acceleration in the Magnetosphere 
and Ionosphere, June 3-7, 1985; Wellesley, Massachusetts. 
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59. N. Singh, H. Thiemann, and R. W. Schunk, Ion beams as seen in simulations of 2- 

D potential structures, presented at the Chapman Conference on Ion Acceleration 
in the Magnetosphere and Ionosphere, June 3-7, 1985; Wellesley, 
Massachusetts. 

60. N. Singh and R. W. Schunk, A possible mechanism for the observed streaming of 

0+ and H + ions at nearly equal speeds in the distant magnetotail, presented at the 
Chapman Conference on Ion Acceleration in the Magnetosphere and Ionosphere, 
June 3-7, 1985; Wellesley, Massachusetts. 

61. R.W. Schunk and J. J. Sojka, A model of the global ionosphere, Invited Review, 

presented at the 1985 North American Radio Science Meeting, June 17-21, 1985; 
Vancouver, British Columbia, Canada. 

62. J. J. Sojka and R. W. Schunk, Theoretical study of the production and decay of 

localized electron density enhancements in the polar ionosphere, presented at the 
1985 North American Radio Science Meeting, June 17-21, 1985; Vancouver, 
British Columbia, Canada. 

63. J. J. Sojka and R. W. Schunk, Theoretical study of anomalously high F-region peak 

altitudes in the polar ionosphere, presented at the 1985 North American Radio 
Science Meeting, June 17-21, 1985; Vancouver, British Columbia, Canada. 

64. C. E. Rasmussen, R. W. Schunk, J. J. Sojka, V. B. Wickwar, O. de la 

Beaujaidiere, J. C. Foster, J. Holt, D. S. Evans, and E. Nielsen, Comparison of 
simultaneous Chatanika and Millstone Hill observations with ionospheric model 
predictions, AGU Spring Meeting, Baltimore, Maryland; EOS, 66, 327, 1985. 

65. N. Singh, A review of mechanisms for the formation of ion conics. Invited Review, 

5th IAGA Scientific Assembly, August 5-17, 1985; Prague, Czechoslovakia. 

66. R. W. Schunk, Models of the terrestrial high latitude ionosphere. Invited Review, 5th 

IAGA Scientific Assembly, August 5-17, 1985; Prague, Czechoslovakia. 

67. A. R. Barakat, R. W. Schunk, T. E. Moore, and J. H. Waite, Ion escape fluxes in 

the terrestrial high latitude ionosphere, 5th IAGA Scientific Assembly, August 5- 
17, 1985; Prague, Czechoslovakia. 

68. A. F. Nagy, A. R. Barakat and R. W. Schunk, Is Jupiter’s ionosphere a significant 

plasma source for its magnetosphere?, 5th IAGA Scientific Assembly, August 5- 
17, 1985; Prague, Czechoslovakia. 

69. H. Thiemann, N. Singh, and R. W. Schunk, Simulation of auroral potential 

structures and related phenomena, 5th IAGA Scientific Assembly, August 5-17, 
1985; Prague, Czechoslovakia. 
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70. R. W. Schunk, Comparison of transport equations for gases and plasmas, and also 

Convergence of the transport equations for gases and plasmas, two one-hour 
lectures presented at the Institut fur Astrophysik und Extraterrestrische 
Forschung, Universitat Bonn, August 15-16, 1985; Bonn, West Germany. 

71. N. Singh, Counterstreaming plasma expansion. Invited Lecture, presented at the 

Nordic Nonlinear Plasma Wave Meeting, Ris0 National Laboratory, August 13- 
16, 1985; Ris0, Denmark. 

72. N. Singh, Simulation of auroral plasma processes, presented at the Department of 

Electron and Plasma Physics, Royal Institute of Technology, August 21, 1985; 
Stockholm, Sweden. 

73. R. W. Schunk, A theoretical global ionospheric F-region model and its predictive 

aspects, Invited Lecture , presented at the SUNDIAL Workshop, September 24- 
26, 1985; McLean, Virginia. 

74. R. W. Schunk and J. J. Sojka, Magnetospheric control of ionospheric processes, 

Invited Review, AGU Fall Meeting, San Francisco, California; EOS, 66, 995, 
1985. 

75. J. J. Sojka and R. W. Schunk, Ionospheric modelling of the GISMOS campaign. 

Invited Talk , AGU Fall Meeting, San Francisco, California; EOS, 66, 994, 1985. 

76. C. E. Rasmussen, R. W. Schunk, and J. J. Sojka, Effects of different convection 

models upon the high-latitude ionosphere, AGU Fall Meeting, San Francisco, 
California; EOS, 66, 996, 1985. 

77. M. D. Bowline, R. W. Schunk, and J. J. Sojka, Theoretical study of the electron 

temperature in the high latitude ionosphere for solar maximum and winter 
conditions, AGU Fall Meeting, San Francisco, California; EOS, 66, 999, 1985. 

78. J. R. Conrad, N. Singh, and R. W. Schunk, EIC waves in a plasma with an ion 

beam and counterstreaming bulk electrons: Waves in the zero-frequency band, 
AGU Fall Meeting, San Francisco, California; EOS, 66, 1034, 1985. 

79. N. Singh, R. W. Schunk, and H. Thiemann, Temporal features of the refilling of a 

plasmaspheric flux tube, AGU Fall Meeting, San Francisco, California; EOS, 66, 
1036, 1985. 

80. H. Thiemann, N. Singh, and R. W. Schunk, Numerical simulations of auroral 

plasma processes, AGU Fall Meeting, San Francisco, California; EOS, 66, 1050, 
1985. 
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81. D. G. Torr, R. W. Schunk, T. L. Killeen, and R. G. Roble, The role of 

thermosphere/ionosphere modeling for the GBOA-CEDAR initiative. Invited 
Talk , AGU Fall Meeting, San Francisco, California; EOS, 66, 990, 1985. 

82. N. Singh and R. W. Schunk, Plasma response to the injection of an electron beam, 

Invited Talk, National Radio Science Meeting, January 13-16, 1986; Boulder, 
Colorado. 

83. R. W. Schunk, Round Table Discussion on Global Ionosphere and Aeronomy Study 

Program, Invited Participant, National Radio Science Meeting, January 13—16, 
1986; Boulder, Colorado. 

84. R. W. Schunk, Space science: High technology now and in the future, Invited 

Lecture, Business and Professional Women’s Organization, April 21, 1986; 
Logan, Utah. 

85. R. W. Schunk, Solar- terrestrial physics: A space age birth, 73 rd Faculty Honor 

Lecture , Utah State University, May 1, 1986; Logan, Utah. 

86. R. W. Schunk, Solar-terrestrial physics, Colloquium, Department of Physics, 

University of Utah, May 8, 1986; Salt Lake City, Utah. 

87. C. E. Rasmussen, P. M. Banks, and R. W. Schunk, Theory of the electrodynamic 

tether, COSPAR, 26 th Plenary Meeting, June 30-July 11, 1986; Toulouse, 
France. 

88. H. Thiemann, R. Grard, N. Singh, and R. W. Schunk, Results of model calculations 

simulating the interaction between the Giotto spacecraft and the impact induced 
secondary plasma, COSPAR, 26 1 * 1 Plenary Meeting, June 30-July 11, 1986; 
Toulouse, France. 

89. J. J. Sojka, R. W. Schunk, and M. D. Bowline, Modelled ionospheric T e profiles at 

high latitudes for possible IRI application, COSPAR, 26 th Plenary Meeting, June 
30-July 11, 1986; Toulouse, France. 

90. R. W. Schunk, An updated theory of the polar wind. Invited Review, COSPAR, 

26 th Plenary Meeting, June 30-July 1 1, 1986; Toulouse, France. 

91. E. P. Szuszczewicz, E. Roelof, R. W. Schunk, B. Fejer, R. Wolf, R. Leitinger, M. 

Abdu, B. M. Reddy, J. Joselyn, P. J. Wilkinson, and R. F. Woodman, A 
cooperative world-wide study of interactive ionospheric processes and their role 
in the transfer of energy and mass in the Sun-Earth system; SCOSTEP 
International Symposium on Solar-Terrestrial Physics, June 30-July 5, 1986; 
Toulouse, France. 
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92. J. J. Sojka, R. W. Schunk, L. A. Frank, J. D. Craven, J. D. Winningham, and J. R. 

Sharber, Ionosphere’s response to an auroral storm based upon the Dynamics 
Explorer SAI and LAPI data bases, SCOSTEP International Symposium on 
Solar-Terrestrial Physics, June 30-July 5, 1986; Toulouse, France. 

93. R. W. Schunk, Interactions between the polar ionosphere and thermosphere, Invited 

Review, SCOSTEP International Symposium on Solar-Terrestrial Physics, June 
30-July 5, 1986; Toulouse, France. 

94. R. W. Schunk, The polar wind. Invited Review, presented at The First Huntsville 

Workshop on Magnetosphere-Ionosphere Plasma Models, October 13-16, 1986; 
Huntsville, Alabama. 

95. A. R. Barakat, R. W. Schunk, T. E. Moore and J. H. Waite, Ion escape fluxes in the 

terrestrial high latitude ionosphere, presented at The First Huntsville Workshop 
on Magnetosphere Ionosphere Plasma Models, October 13-16, 1986; Huntsville, 
Alabama. 

96. R. W. Schunk, Ionosphere-thermosphere modelling: Current status. Invited Review, 

presented at the International Symposium on Large-Scale Processes in the 
Ionospheric-Thermospheric System, December 2-5, 1986; Boulder, Colorado. 

97. C. E. Rasmussen, R. W. Schunk, J. J. Sojka, V. B. Wickwar, O. de la 

Beaujardiere, J. C. Foster, and J. M. Holt, Comparison of simultaneous 
Chatanika and Millstone Hill temperature measurements with ionospheric model 
predictions, presented at the International Symposium on Large-Scale Processes 
in the Ionospheric-Thermospheric System, December 2-5, 1986; Boulder, 
Colorado. . 

98. J. J. Sojka and R. W. Schunk, Asymmetries in the plasma characteristics of the 

conjugate high-latitude ionospheres, presented at the International Symposium on 
Large-Scale Processes in the Ionospheric-Thermospheric System, December 2-5, 
1986; Boulder, Colorado. 

99. C. E. Rasmussen and R. W. Schunk, An E and F region density model for obtaining 

ionospheric conductivity, presented at the International Symposium on Large- 
Scale. Processes in the Ionospheric-Thermospheric System, December 2-5, 1986; 
Boulder, Colorado. 

100. E. P. Szuszczewicz, E. C. Roelof, R. W. Schunk, B. G. Fejer, R. A. Wolf, R. 

Leitinger, M. A. Abdu, B. M. Reddy, J. C. Joselyn, P. J. Wilkinson, and R. F. 
Woodman, SUNDIAL: A continuing program of global-scale modeling and 
measurement of ionospheric responses to solar, magnetospheric, and 
thermospheric controls, presented at the International Symposium on Large-Scale 
Processes in the Ionospheric-Thermospheric System, December 2-5, 1986; 
Boulder, Colorado. 
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101. R. W. Schunk, Ionospheric physics: The impact of plasma transport. Invited Talk, 

AGU Fall Meeting, San Francisco, California; EOS, 67, 1 126, 1986. 

102. H. G. Demars and R. W. Schunk, Comparison of solutions to bi-Maxwellian and 

Maxwellian transport equations for subsonic flows, AGU Fall Meeting, San 
Francisco, California; EOS, 67, 1136, 1986. 

103. A. R. Barakat and R. W. Schunk, Stability of the polar wind, AGU Fall Meeting, 

San Francisco, California; EOS, 67, 1 136, 1986. 

104. /. J. Sojka and R. W. Schunk, Theoretical study of the high latitude ionosphere’s 

response to multi-cell convection patterns, AGU Fall Meeting, San Francisco, 
California; EOS, 67, 1137, 1986. 

105. R. J. Sica, C. E. Rasmussen, and R. W. Schunk, Can the high latitude ionosphere 

support large field-aligned ion drifts?, AGU Fall Meeting, San Francisco, 
California; EOS, 67, 1137, 1986. 

106. R. W. Schunk and J. J. Sojka, The lifetime and transport of severe ionospheric 

disturbances, AGU Fall Meeting, San Francisco, California; EOS, 67, 1137, 
1986. 

107. C. E. Rasmussen and R. W. Schunk, Ionospheric convection driven by NBZ 

currents, AGU Fall Meeting, San Francisco, California; EOS, 67, 1162, 1986. 

108. C. E. Rasmussen and R. W. Schunk, Ionospheric convection driven by asymmetries 

in NBZ currents, National Radio Science Meeting, January 12-15, 1987; 
Boulder, Colorado. 

109. R. W. Schunk, Implications of the SUNDIAL data on first-principles global-scale 

models, presented at the Third SUNDIAL Workshop, February 24-27, 1987; La 
Jolla, California. 

110. R. W. Schunk and J. J. Sojka, Ionospheric features induced by magnetospheric 

processes, Invited Talk, presented at the International Symposium on Quantitative 
Modeling of Magnetosphere-Ionosphere Coupling Processes, March 9-13, 1987; 
Kyoto, Japan. 

111. E. P. Szuszczewicz, E. Roelof, R. W. Schunk, B. Fejer, R. Wolf, M. Abdu, J. 

Joselyn, B. M. Reddy, P. Wilkinson, R. Woodman, and R. Leitinger, 
SUNDIAL: A global study of ionospheric processes and their roles in the 
transfer of energy and mass in the sun-earth system, AGU Spring Meeting, 
Baltimore, Maryland; EOS, 68, 364, 1987. 
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1 12. R. W. Schunk and E. P. Szuszczewicz, First-principle and empirical modeling of the 

global-scale ionosphere. Invited Talk , AGU Spring Meeting, Baltimore, 
Maryland; EOS, 68, 364, 1987. 

113. P. J. Wilkinson, R. W. Schunk, R. Hanbaba, and E. P. Szuszczewicz, 

Interhemispheric comparison of SUNDIAL F-region data with global-scale 
ionospheric models, AGU Spring Meeting, Baltimore, Maryland; EOS, 68, 364, 
1987. 

114. M. A. Abdu, B. M. Reddy, G. O. Walker, R. Hanbaba, J. H. A. Sobral, B. G. 

Fejer, R. W. Schunk, R. F. Woodman, and E. P. Szuszczewicz, Processes in 
the quiet and disturbed equatorial low latitude ionosphere: SUNDIAL campaign 
1984, AGU Spring Meeting, Baltimore, Maryland; EOS, 68, 364, 1987. 

115. R. W. Schunk, Polar wind theory. Invited Review, presented at the 1987 Cambridge 

Workshop in Theoretical Geoplasma Physics on “Ionosphere-Magnetosphere- 
Solar Wind Coupling Processes,” July 28-August 1, 1987; Cambridge, 
Massachusetts. 

116. R. W. Schunk and E. P. Szuszczewicz, First-principal and empirical modelling of the 

global- scale ionosphere, presented at the International Union of Geodesy and 
Geophysics, August 9-22, 1987; Vancouver, British Columbia, Canada. 

117. R. J. Sica, C. E. Rasmussen, and R. W. Schunk, Can the high-latitude ionosphere 

support large field-aligned ion drifts?, presented at the International Union of 
Geodesy and Geophysics, August 9-22, 1987; Vancouver, British Columbia, 
Canada. 

118. C. E. Rasmussen and R. W. Schunk, Ionospheric convection driven by NBZ 

currents, presented at the International Union of Geodesy and Geophysics, 
August 9-22, 1987; Vancouver, British Columbia, Canada. 

119. A. R. Barakat and R. W. Schunk, Stability of the polar wind: Linear theory and 

simulation, presented at the International Union of Geodesy and Geophysics, 
August 9-22, 1987; Vancouver, British Columbia, Canada. 

120. C. E. Rasmussen and R. W. Schunk, An E- and F-region density model for 

obtaining ionospheric conductivity, presented at the International Union of 
Geodesy and Geophysics, August 9-22, 1987; Vancouver, British Columbia, 
Canada. 

121. C. E. Rasmussen, J. J. Sojka, R. W. Schunk, V. B. Wickwar, O. de la 

Beaujardiere, J. Foster, and J. Holt, Comparison of simultaneous Chatanika and 
Millstone Hill temperature measurements with ionospheric model predictions, 
presented at the International Union of Geodesy and Geophysics, August 9-22, 
1987; Vancouver, British Columbia, Canada. 
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122. J. J. Sojka, New results in modelling the global and high latitude ionospheric 

dynamics, presented at the International Union of Geodesy and Geophysics, 
August 9-22, 1987; Vancouver, British Columbia, Canada. 

123. R. W. Schunk, Ionosphere-thermosphere interaction modelling; Current status, 

Invited Paper , AGU Fall Meeting, San Francisco, California; EOS, 68, 1382, 
1987. 

124. H. G. Demars and R. W. Schunk, Temperature anisotropies in the ionosphere, AGU 

Fall Meeting, San Francisco, California; EOS, 68, 1382, 1987. 

125. J. J. Sojka and R. W. Schunk, A model study of how electric field structures affect 

the polar cap F-region, AGU Fall Meeting, San Francisco, California; EOS, 68, 
1385, 1987. 

126. A. R. Barakat and R. W. Schunk, Effect of cusp H + and 0 + beams on the stability of 

the polar wind, AGU Fall Meeting, San Francisco, California; EOS, 68, 1386, 

1987. 

127. W.-H. Yang and R. W. Schunk, A model of solar wind streams, AGU Fall Meeting, 

San Francisco, California; EOS, 68, 1410, 1987. 

128. C. E. Rasmussen and R. W. Schunk, Multi-fluid approach to study the refilling of 

plasmaspheric flux tubes, AGU Fall Meeting, San Francisco, California; EOS, 
68, 1421, 1987. 

129. J. Liu, C. E. Rasmussen, R. J. Sica, R. W. Schunk, and A. D. Richmond, Motion 

of ionospheric flux tubes during a period of substorm activity, AGU Fall 
Meeting, San Francisco, California; EOS, 68, 1430, 1987. 

130. R. W. Schunk and J. J. Sojka, Modelling ionospheric density structures. Invited 

Talk, presented at the National Radio Science Meeting, January 5-8, 1988; 
Boulder, Colorado. 

131. R. W. Schunk, A. R. Barakat, H. G. Demars, T.-Z. Ma, C. E. Rasmussen, J. J. 

Sojka, H. Thiemann, and W.-H. Yang, The flow of plasma in the solar-terrestrial 
environment. Invited, AGU Spring Meeting, Baltimore, Maryland; EOS, 69, 
423, 1988. 

132. C. E. Rasmussen and R. W. Schunk, A three-dimensional time-dependent model of 

the plasmasphere, AGU Spring Meeting, Baltimore, Maryland; EOS, 69, 428, 

1988. 

133. R. W. Schunk and J. J. Sojka, A three-dimensional time-dependent model of the 

polar wind, AGU Spring Meeting, Baltimore, Maryland; EOS, 69, 429, 1988. 
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134. A. R. Barakat and R. W. Schunk, Stability of energetic ion outflows from the 

ionosphere, AGU Spring Meeting, Baltimore, Maryland; EOS, 69, 429, 1988. 

135. J. J. Sojka and R. W. Schunk, Modelling ionospheric density structures and dynamic 

auroral features, AGU Spring Meeting, Baltimore, Maryland; EOS, 69, 431, 
1988. 

136. T.-Z. Ma and R. W. Schunk, A two-dimensional model of plasma expansion in the 

ionosphere, AGU Spring Meeting, Baltimore, Maryland; EOS, 69, 431, 1988. 

137. R. J. Sica, C. E. Rasmussen, and R. W. Schunk, Maximum field-aligned plasma 

velocities in the f-region, AGU Spring Meeting, Baltimore, Maryland; EOS, 69, 
443, 1988. 

138. W.-H. Yang and R. W. Schunk, Modelling solar wind streams and the neutral 

current sheet, AGU Spring Meeting, Baltimore, Maryland; EOS, 69, 436, 1988. 

139. R. W. Schunk and J. J. Sojka, Modelling ionospheric density structures, AGARD 

Symposium on “Ionospheric Structure and Variability on a Global Scale and 
Interaction with Atmosphere, Magnetosphere,” May 16-20, 1988; Munich, West 
Germany. 

140. R. W. Schunk, Polar wind tutorial. Invited Lecture, presented at the 1988 Cambridge 

Workshop in Theoretical Geoplasma Physics on “Polar Cap Dynamics and High 
Latitude Ionospheric Turbulence,” June 13-17, 1988; Cambridge, 
Massachusetts. 

141. R.W. Schunk, Magnetosphere-ionosphere-thermosphere coupling processes. Invited 

Talk, presented at the SCOSTEP special session on “Solar-Terrestrial Energy 
Program (STEP): Major Scientific Problems,” July 18-29, 1988; Helsinki, 
Finland. 

142. R. W. Schunk, Response of the ionosphere-thermosphere system to magnetospheric 

forcing, Invited Talk, COSPAR 27 th Plenary Meeting, July 18-29, 1988; 
Helsinki, Finland. 

143. A. R. Barakat and R. W. Schunk, Outflow of low energy ions in the terrestrial high- 

latitude ionosphere, Invited Talk, COSPAR 27 th Plenary Meeting, July 18-29, 
1988; Helsinki, Finland. 

144. A. R. Barakat and R. W. Schunk, Stability of the polar wind in the presence of cusp 

H+ and 0+ ions, COSPAR 27 th Plenary Meeting, July 18-29, 1988; Helsinki, 
Finland. 
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145. H. G. Demars and R. W. Schunk, Solutions to bi-Maxwellian transport equations for 

the solar wind, American Physical Society Topical Conference on Plasma 
Astrophysics, September 19-23, 1988; Santa Fe, New Mexico. 

146. R. W. Schunk and W.-H. Yang, Modeling high-speed solar wind streams, American 

Physical Society Topical Conference on Plasma Astrophysics, September 19-23, 
1988; Santa Fe, New Mexico. 

147. W.-H. Yang, Force-free magnetic field model of accretion disk, American Physical 

Society Tropical Conference on Plasma Astrophysics, September 19-23, 1988; 
Sante Fe, New Mexico. 

148. J. A. Klimchuk, P. A. Stuirock, and W.-H. Yang, Coronal magnetic fields produced 

by photospheric shear, 172 nd American Astronomical Society, Kansas City, 
1988. 

149. R. W. Schunk and J. J. Sojka, Global polar wind variations during changing 

magnetospheric conditions, AGU Fall Meeting, San Francisco, California; EOS, 
69, 1340, 1988. 

150. J. J. Sojka and R. W. Schunk, Theoretical study of the seasonal behavior of the 

global ionosphere at solar maximum, AGU Fall Meeting, San Francisco, 
California; EOS, 69, 1349, 1988. 

151. H.G. Demars and R. W. Schunk, Solutions to bi-Maxwellian transport equations for 

the solar wind, AGU Fall Meeting, San Francisco, California; EOS, 69, 1364, 
1988. 

152. W.-H. Yang and R. W. Schunk, Modelling high-speed solar wind streams, AGU 

Fall Meeting, San Francisco, California; EOS, 69, 1358, 1988. 

153. R. J. Sica, R. W. Schunk, and P. Wilkinson, An empirical study of the undisturbed 

mid-latitude ionosphere using simultaneous, multiple site ionosonde 
measurements, AGU Fall Meeting, San Francisco, California; EOS, 69, 1350, 
1988. 

154. C. E. Rasmussen and R. W. Schunk, Comparison of plasmaspheric measurements 

with a three-dimensional time-dependent model of the plasmasphere, AGU Fall 
Meeting, San Francisco, California; EOS, 69, 1347, 1988. 

155. R. W. Schunk, Modelling highly non-equilibrium space plasmas. Colloquium at the 

National Center for Atmospheric Research, January 5, 1989; Boulder, Colorado. 

156. R. W. Schunk, Current status of numerical ionospheric modelling. Invited Talk, 

AGU Spring Meeting, Baltimore, Maryland, EOS, 70, 406, 1989. 
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157. V. B. Wickwar and R. W. Schunk, The high-latitude F region: Observation and 

theoiy, Invited Talk, presented at the Fourth EISCAT Scientific Workshop, June 
5-9, 1989; Sigtuna, Sweden. 

158. J. J. Sojka, Ionospheric physics, Invited Talk, Summer Institute on Atmospheric 

Sciences, sponsored by NASA Goddard Space Flight Center, June 21, 1989, 
NASA, Goddard, Greenbelt, Maryland. 

159. T.-Z. Ma and R. W. Schunk, A 3-D model for plasma clouds in the ionosphere, 

IAGA Scientific Assembly, July 24-August 4, Exeter, England; I AG A Bulletin, 
53, 393, 1989. 

160. R. W. Schunk and E. P. Szuszczewicz, Plasma expansion phenomena: Comparison 

of characteristics predicted by small-scale and macroscopic formulations, IAGA 
Scientific Assembly, July 24-August 4, Exeter, England; IAGA Bulletin, 53, 
393, 1989. 

161. R. W. Schunk and /. J. Sojka, Global polar wind variations driven by 

magnetospheric processes, IAGA Scientific Assembly, July 24-August 4, 
Exeter, England; IAGA Bulletin, 53, 309, 1989. 



